Field-induced magnetoelastic instabilities in antiferromagnetic molecular wheels 
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The magnetic torque of the antiferromagnetic molecular wheel CsFeg was studied down to 50 mK 
and in fields up to 28 T. Below ca. 0.5 K phase transitions were observed at the field-induced 
level- crossings (LCs). Intermolecular magnetic interactions are very weak excluding an explanation 
in terms of field-induced magnetic ordering. A magneto-elastic coupling was considered. A generic 
model shows that the wheel structure is unconditionally unstable at the LCs, and the predicted 
torque curves explain the essential features of the data well. 

PACS numbers: 33.15.Kr, 71.70.-d, 75.10.Jm 
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Antiferromagnetic (AF) molecular wheels attracted 
huge attention recently because of their peculiar quantum 
properties P, S S H ■ These molecules are character- 
ized by a ring-like arrangement of magnetic metal ions; 
the ferric wheel [CsFe8{N(CH2CH20)3}8]Cl, or CsFeg 
0, studied in this work is shown in Fig.^a). The ions 
within a wheel experience AF nearest-neighbor Heisen- 
berg interactions, and the molecule's ground state at zero 
magnetic field is nonmagnetic with total spin 5* = 0. The 
next higher lying states belong to S" = 1,2, etc. In a 
magnetic field these states split due to the Zeeman in- 
teraction, leading to a series of level crossings (LCs) at 
characteristic fields at which the ground state changes 
from the S* = 0, M = level to the S" = 1, M = -1 
level, the S = 2, M = —2 level, and so on, see inset 
of Fig. ^b) P, 0. Hi- The magnetization curve at low 
temperatures thus exhibits a staircase-like field depen- 
dence, with a step at each LC. In this work we report 
field-dependent measurements of the magnetic torque on 
CsFeg, which show clear indications of phase transitions 
at the LCs at low temperatures. 

The observed anomalies at the LCs could be due to 
weak magnetic interactions between the molecules in the 
sample (mier molecular interactions). Since at low tem- 
peratures AF wheels behave like dimers 3, this would 
place CsFeg in the context of weakly-interacting dimer 
compounds, such as TlCuCls, which may exhibit Bose- 
Einstein condensation of magnons 

'91,0. 

Observation of 

such phenomena in a crystal of molecular wheels would 
be of great interest, but in CsFeg intermolecular interac- 
tions are very weak excluding such a scenario. 

The degeneracy at the LCs suggests a Spin-Peierls type 
of effect as an alternative, in which the degeneracy is 
lifted due to a coupling of the spin system to the lattice. 
CsFeg is well described by the spin Hamiltonian 0, 



^S, -S.+i + Sg-Si -K^MBS-B-^i/A, (1) 



magnetic anisotropy, mostly due to ligand-field and dipo- 
lar interactions (the magnetic anisotropy is well described 
by a term i'X^i^i ^1 z with D — -0.56 K and the wheel 
axis z). Si is the spin operator of the ?th ion with spin 
s = 5/2. Recent theoretical studies on isotropic AF spin 
rings {Ha — 0) yielded a magneto-elastic (ME) instabil- 
ity in zero field for sufficiently small elastic constants, but 
only for rings with s = 1/2 [l^ Il3|. These conclusions 
hold also with applied magnetic fields since a modulation 
of the exchange constants along the ring cannot lift the 
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which includes the Heisenberg and Zeeman terms {J — - 
20.6 K, g = 2), and a term Ha describing a weak uniaxial 



FIG. 1: a) Torque vs. field, t{B), for a single crystal of 2 at 
55, 100, 200, 300, and 400 mK (online colors: black, red, blue, 
green, magenta) (95 = -3.3°). The right inset shows dr/dB 
at the first LC, arrows indicate the onset of anomalies. The 
left inset shows the structure of CsFeg. b) r(_B) for a single 
crystal of 1 at 60 and 700 mK (online colors: black, red) {ip 
— 93.6°). The right inset shows the data as dr/dB. The left 
inset shows a calculated energy spectrum as function of field 
neglecting magnetic anisotropy (J = —20.6 K, arrows mark 
the first two LCs, the involved levels are labeled by S, M). 



2 



degeneracy at the field-induced LCs. However, as will 
be shown below, the situation changes markedly in the 
presence of a magnetic anisotropy, as in CsFcg. A simple, 
generic model is introduced, which demonstrates that at 
the LCs the ring structure is unstable against distortions 
for any value of the spin-phonon coupling. The model 
reproduces the characteristic features of the data well, 
suggesting a ME origin of the anomalies at the LCs in 
CsFeg. 

Single crystals of CsFeg were prepared as in Rcf. 0, 
but crystallized either from a mixture of CHCI3 
and CH2CI2 by pentane vapor diffusion yielding 
CsFeg-SCHCla-O.SHaO (1) [H, or from ethanol by 
diethyl-ether vapor diffusion yielding CsFe8-8C2H50H 
(2) 1 (2) crystalhzes in the space group P21/n 

(P4/n) and the molecules exhibit approximate (crystal- 
lographic) C4 symmetry. The magnetic parameters, J 
and £>, of the molecules are not affected by the solvents 
but m<ermolecular exchange interactions, if present, 
should be strongly affected. The magnetic torque r was 
measured with a CuBe cantilever inserted into the MIO 
magnet at the Grenoble High Magnetic Field Laboratory 
equipped with an Oxford 3He/4He dilution fridge. Back- 
ground signals were negligible; only raw data are shown 
here. In total six single-crystal samples were investigated. 

Figure n presents the field dependence of the torque t 
of CsFeg at various temperatures for fields close to the 
uniaxial axis and almost perpendicular to it {(p denotes 
the angle between field and uniaxial axis z). At the 
highest temperatures, the curves exhibit the expected 
behavior: At the LC fields the torque curves display 
steps broadened by the effect of, e.g., temperature (a 
plot dr/dB vs. B shows Gaussian- like peaks with widths 
r). For CsFeg, the first three LCs are observed in fields 
up to 28 T. Normally, the steps become sharper with 
lower temperatures, corresponding to decreasing widths 
r, but in CsFeg, in contrast, a very different behavior 
is observed at low temperatures. For fields close to the 
uniaxial axis, « 0°, a dome-shaped contribution to the 
torque, centered at the LC fields, appears. This is ap- 
parent in Fig.IHa) for the LCs at 18 T and 25.5 T, but 
also the first LC at 11 T shows the anomaly, as seen in 
dr/dB [inset of Fig.^a)]. Anomalies were also observed 
for close to perpendicular fields, ip w 90°, Fig.Q^b). Here 
the torque exhibits a linear field dependence at the first 
LC at 7.7 T, but also at the higher LCs at 16.5 T and 
24 T anomalous behavior is evident from the inset of 

Fig.inb). 

Figure ISfa) presents dr/dB near the first LC for ip k, 
90° as determined from measurements at several tem- 
peratures from 55 mK to 1 K. At the higher tempera- 
tures, the curves exhibit the usual Gaussian-like field de- 
pendence, but below a critical temperature Tc, of about 
0.65 K for the shown sample, a plateau corresponding to 
a linear field dependence in r{B) emerges. The temper- 
ature dependencies of the lower and upper critical fields. 
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FIG. 2: a) dr/dB at the first LC for a single crystal of 
1 for several temperatures from 55 mK to 1 K ((^ ~ 93°). 
b) B-T phase diagram as derived from the data shown in 
panel a) (lines are guides to the eyes). The closed symbols 
indicate the critical fields, the open symbols the half-maximun 
fields for T > Tc. The dashed line indicates the field of the 
S' = — > 5 = 1 LC. The dash-dotted lines indicate the half- 
maximum fields as expected for thermal broadening. 

as defined by the kinks in dr/dB, are plotted in Fig.E^b) 
(the difference will be denoted as ABq). The figure also 
displays the maximum and half-maximum field values of 
the Gaussian- like curves for T > Tc- 

All investigated crystals showed the anomalies at the 
LCs, with a variation in Tc (0.35 to 0.65 K) and Ai?o 
(1.5 to 3 T). Smaller ABq corresponded to smaller 
and vice versa. One origin for the variation seems to 
be a dependence of Tj, and ASq on the field orientation 
indicated by the data. Within experimental resolution, 
a hysteresis was not detected. 

The above torque data clearly demonstrates a phase 
transition at the LCs in CsFeg. The similarity of the 
torque curves in Fig.^b) with magnetization curves ob- 
served for systems exhibiting field-induced magnetic or- 
der ^3 suggests to assign the anomalies in CsFeg to 
intermolecular magnetic interactions. However, this is 
an unlikely scenario because i) intermolecular magnetic 
dipole-dipole interactions are very weak (even optimistic 
estimates yield values < 10 mK), ii) in the crystal struc- 
tures of 1 and 2 the individual molecules are well sep- 
arated and no exchange pathways exist, iii) exchange 
interactions would strongly depend on the solvent and 
the details of the crystal packing, in contrast to the 
observation of similar critical field ranges and temper- 
atures in crystals of 1 and 2, and iv) no indications 
of intermolecular interactions were found in the compa- 
rable systems [NaFe6{N(CH2CH203}6]Cl-6CHCl3 down 
to 0.3 K, [Fe6{N(CH2CH203}6]-6MeOH down to 0.2 K, 
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andJNaFfefi(OCH3)i2(Ci7Hi504)6]C104 down to 40 mK 
[l^ I20I El| . Intermolecular interactions on the order of 
several 100 mK are not apparent in CsFcg. 

In view of the degeneracy at the LCs a ME effect might 
be suggested. For isotropic AF spin ring s such an effect 
is unrealistic, at least for spins > 1/2 |12LIi^. but, as will 
be shown below, with an additional magnetic anisotropy 
the ring structure becomes unconditionally unstable at 
the LCs: For a distorted ring the magnetic anisotropy 
allows for a mixing of the levels at a LC, resulting in 
an avoided LC with a gap A and hence a lowering of the 
ground-state energy by A/2, which is proportional to the 
modulation of the anisotropy constants and thus in first 
approximation linear in the distortion. 

Following Ref. [l^l, the low-T magnetic behavior near 
a LC is well described by the two-level Hamiltonian 

f es A/2 
I A/2 65+1 



s+i 



(2) 



where S and S + 1 index the two levels |5, — S") and \S + 
1, —S — 1) involved in the level crossing. Here, a rotated 
axis system with the quantization axis parallel to the field 
B is used, and 15*,^/) refers to the labels of the total 

spin operator in the rotated frame, S and S'^. es{B, ip) 
describes the field and angle dependence of the levels 
without a mixing. In first order, one obtains es{B,Lp) — 
—bS + Asif) with the reduced field b — gfJ-sB, where 
As{(p) accounts for the exchange interactions and the 
zero-field splittings produced by Ha 0j Ull' The LC 
field is given by bo{f) = A5+1 — Ag. A level mixing 
at the LC is included in the model via A (1^9), which in 
first order is given by A/2 = {S, -S\H'j^\S +1,-S - 1) 
and is thus independent of the magnetic field (iJ^ is Ha 
expressed in the rotated reference frame). 

The symmetry of a non-distorted ring (Cg symmetry 
of the spin Hamiltonian) prohibits a mixing of the levels, 
hence A = 0, ^ structural distortion induces 

a modulation of the exchange and anisotropy constants, 
which affects both A5 and A [the ME coupling due to 
A5 (A) will be called diagonal (non-diagonal)]. Ag is not 
affected in first order by these modulations as it is only 
sensitive to the average of the exchange and anisotropy 
constants ,11]. Thus, A5 varies as A5 cc x^, where x is a 
parameter describing the structural distortion (as usual it 
has been assumed that the modulation of the parameters 
is linear in the distortion). This is in accordance with 
Refs. 0,^3, and leads to similar conclusions. For the 
non-diagonal ME coupling, however, since it is sensitive 
to the amplitude of the modulations, one finds A oc a; 
[we write A(x) = ax with the ME coupling constant 
a]. The gain in magnetic energy is now linear in the 
distortion resulting in an unconditional ME instability 
|23| . In this model, the microscopic details of the non- 
diagonal ME coupling are lumped into the parameter a. 
Questions concerning, e.g., the relevant distortion mode 
thus remain unanswered. 




FIG. 3: Field dependence of a) the order parameter Ao, b) 
the torque contribution 5ti, and c) the torque contribution 
St2 around a LC at bo [the plots were normalized to o? /{2k), 
dbo/dtp, and —dAo/dif, respectively]. 



The change in ground-state energy due to an opening 
of a gap at the LC is easily calculated. Including the 
elastic energy, the potential V{x) of the total system is 



Vix) 



'-V{b-bor + A{xr + -\b- 



bo\ + -kx^, (3) 



where k is the spring constant. In the adiabatic approxi- 
mation, which neglects the kinetic energy of the phonons 
(the validity of this approach is discussed in Ref. [l^), 
the equilibrium distortion xo is given by the minimum of 
V{x), yielding the condition 



(b-bo)' 



(4) 



with Ao = A(a;o). Accordingly, as a function of magnetic 
field, the order parameter Aq describes a semi circle of 
radius a'^/{2k) around fopi see Fig. (SJa), and the system 
exhibits a spontaneous distortion for fields in between 
60 ± be, with be = a^/{2k). 

The model permits an analytical calculation of the 
torque profiles at low temperature. Inserting eq. Q into 
eq. (4) of Ref. jl^ , the change of the torque 6t for mag- 
netic fields in the critical region 60 i be is obtained as 



1 



(b-bo) 



1 dAp Aq 

2 dip be 



(5) 



The torque signal consists of two contributions. The first 
term, 6ti , describes a linear increase of the torque from 
zero at 6 < 60 — to dbo/dip at b > bo + be [Fig. Elb)], 
while the second term, St2 , is proportional to the square 
of the order parameter, St2 oc Aq [Fig. I^c)]. The rele- 
vance with the experimental data is obvious: the dome- 
shaped torque contribution produced by St2 resembles 
the behavior at the LCs for w 0° , Fig. da) , and the 
slope-like contribution of 5ti that for ip « 90° , Fig. GJb) . 

In general the torque signal in the critical region is 
a combination of a slope-like and a dome-shaped curve, 
with angle dependent relative weights. The dependence 
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of the LC field on turns out to be bo{ip) oc 3 cos^ <p— 1 
llH , so that the contribution of the slope- hke part to the 
torque, which is controlled by the factor dbo/dip, varies 
as cos sin (/J. The variation of the contribution of the 
dome-shaped part with angle is more complex, since the 
various magnetic anisotropy terms of possible relevance 
may result in very different angle dependencies of Aq . An 
anisotropy term Hjj = J2 ^A^iz ~ Si{Si + l)/3], for in- 
stance, gives rise to A oc sin(c^) cos{if). A Dzyaloshinsky- 
Moriya interaction Hjji^i = ^ x S^+i), on the other 
hand, which is likely to arise in the course of a structural 
distortion because of the lowered symmetry of the ring 
[l8| . varies as sin((/3), so that dAo/d(p is important at 
small angles but negligible near 90°. This could explain 
a dominance of the dome-shaped contribution for parallel 
fields, and of the slope-like contribution for perpendicular 
fields, as observed. The details are not yet understood, 
and more studies are clearly needed. However, the pro- 
posed scenario is capable of explaining the different find- 
ings for nearly parallel and perpendicular fields, which 
in our opinion supports the idea of a ME origin of the 
observed anomalies in CsFcs. 

With fc = 10 N/m (corresponding to a typical phonon 
frequency of = 100 cm^^ and a reduced mass /i = 
14), a be of about 1 T implies the reasonable value a 
= 3 meV/A (in spin-Peierls systems, e.g., the coupling 
constant is ~ 10\J\/do with do ~ 3.5 A). Concerning 
the stability of the distortions against thermal fluctua- 
tions, the critical temperature may be estimated from 
l/2{d^V/dx'^)xl « l/2kBTe, i.e., keTc « 6c/2. A crit- 
ical field of about 1 T then suggests Tc ~ 1 K, which 
is on the order of the observed values. Here however it 
should be noted that for systems, which exhibit uncondi- 
tional lattice instabilities, BCS type of relations between 
order parameter and are found 0. Also, strain ef- 
fects between the molecules in the crystal might result in 
cooperativity which would help to stabilize the distorted 
phase ^2A\ . A realistic theory thus might have to include 
not only the optic but also the acoustic phonons. 

In conclusion, we have studied the field dependence of 
the magnetic torque for the AF molecular wheel CsFeg 
and observed anomalies at the level-crossing fields at low 
temperatures. With respect to their explanation, sev- 
eral mechanisms were considered. Magnetic interactions 
between different molecules are very weak and are thus 
unlikely to cause the anomalies. As a second possibil- 
ity, magneto-elastic instabilities were discussed. Indeed, 
by introducing a generic model, it has been shown that 
a non-diagonal magneto-elastic coupling due to a mag- 
netic anisotropy induces structural instabilities at the 
LCs. The predicted torque curves allowed to explain 
the generic features of the experimental data, in partic- 
ular the unusual dome-shaped parts for magnetic fields 
close to the uniaxial axis. The current work necessarily 
could not answer all questions, and more experimental 
as well as theoretical work is needed. Of most impor- 



tance would be structural measurements. On the one 
hand, they would allow a direct test of the above sce- 
nario, and on the other hand, would yield information 
about the relevant distortion modes - a crucial input for 
any microscopic model to be developed. 

For a number of ferric wheels, evidence for gaps at the 
LCs has been reported 0, E, 0, |^. In these works, 
thermodynamic data were found to be better explained 
by assuming avoided LCs (with temperature-independent 
gaps). In this context it is interesting to note that the 
anomalies observed in CsFeg below Tc are announced by 
an excessive broadening of the torque steps above Tc, see 
Fig. [S^b). This mimics avoided LCs, and it will be thus 
interesting to see whether the earlier reports of avoided 
LCs were not in fact first hints of the above anomalies. 
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